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1 Introduction 



Studies of the structure of atoms and ions (ultracold to relativistic included) as well as their two- 
body interactions with photons, electrons and other particles require accurate methods for the 
description of such objects. 

In order to obtain accurate values of atomic quantities it is necessary to account for relativistic 
Q and correlation effects. Relativistic effects may be taken into account as Breit-Pauli corrections 
or in a fully relativistic approach. In both cases for complex atoms and ions, a considerable part 
J>^, of the effort must be devoted to integrations over spin-angular variables, occurring in the matrix 
' elements of the operators under consideration. 

Many existing codes for integrating are based on a scheme by Fano [1]. The integrations over 
^ . spin-angular variables in this case constitute a considerable part of the problem, especially when 
<^ atoms with many open shells are treated, and the operators are not trivial. In the papers of 
Gaigalas et al [2, 3], the efficient approach for finding matrix elements of any one- and two-particle 
atomic operator between complex configurations is suggested. It is free of shortcomings of previous 
approaches. This approach allows one to generate fairly accurate databases of atomic parameters 
(Froese Fischer et al [4, 5]). 

Further development of the approaches by Gaigalas et al [2, 3] for the spin-spin and spin-other- 
orbit relativistic corrections in the Breit-Pauli approximation is presented in this poster. 



o 

O 

^ 2 Matrix Elements Between Complex Configurations 

According to the approach of Gaigalas et al [3], a general expression of the submatrix element for 
any two-particle operator between functions with u open shells, can be written as follows: 
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where A = {LiSi,LjSj,L[S[,L'-S'A is the array of bra function shells' terms, and similarly 



i) 

for A*^^* and X = ls. 
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Thus, to calculate the spin-angular part of a submatrix element of this type, one has to obtain: 

1. The recoupling matrix i? ^Aj, Aj, A^, A^, A*""", A'^^*, , which has an analytical expression in 
term of just 6j- and 9j-coefficients. 

2. Submatrix elements T ^njAj, rijAj, n^A^, n^A^, A'"'", A*^'^*, H, , for tensorial products of cre- 
ation/annihilation operators that act upon a particular electron shell. So, all the advantages 
of tensorial algebra and quasispin formalism (Rudzikas [6]) may be efficiently exploited in the 
process of their calculation. 

3. Phase factor A, 

4. 0' (riiXi, nj\j,n^X[, n'jXpE^ , which is proportional to the two-electron submatrix element of 
operator G. 

Further development of this approach for the spin-spin and spin-other-orbit relativistic correc- 
tions in the Breit-Pauli approximation is presented in the following section. 

3 The Spin-Spin and Spin-Other-Orbit Operators 

The spin-spin operator iJ** itself contains tensorial structure of two different types, summed over 
k: 



j^ss — r^(fc+lfe-12,112) I ^(A;-lfe+12,112)' 

/ ss ss 

k 



(2) 



Their submatrix elements are: 



{mXiUjXj \H^l+^^-^'^^^^'^'^\ni,Ximj,Xj>j = -^^J {2k + 2.)^^'^ x 



(3) 



(niXiUjXj ||//i^-i'=+^2,ii2)||^.,_^.,^^.,^^.,j ^ A^^(^2k + 3f^ X 
X [li \\C^''~^ k') (ij ||C('=+^)|| If') N''-^ {njljnili,nflfnidi>) , 



(4) 



where we use a shorthand notation (2k + 3)^^^ = {2k + 3) {2k + 2) {2k + 1) {2k) {2k - 1) and radial 
integral in (3), (4) is defined as in Glass and Hibbert [7]: 



a 



'0^0 I I 

where e{x) is a Heaviside step-function. 



A^'^ {uiliUjlj^niiliinjilji) = 

k 

Pi {ri) Pj (r2) -^e{ri - r2)Pi' {n) Pf (rs) dndr2,, 



(5) 



2 



1; for X > 0, 
0; for x<0. 



(6) 



The spin- other- orbit operator H^°° itself contains tensorial structure of six different types, 
summed over k: 



TTsoo — \^ ro-(A;-lfel,101) I r7-(fe-lfcl,011) I o-(fefel,101) , r7-(feA;l,011) , o-(fe+lfel,101) i r7-(fe+lfel,011) 

— / J SOO ' SOO ' SOO ' SOO ' SOO ' SOO 



Their submatrix elements are: 



(7) 



(niXiUjXj \\h^oo^'''^'''"'^'^^\\ ni'Xi'UfXj,^ = 2-2''^ {{2k - 1) (2A; + 1) x 
X {li -k + l){k-li + kf) {k + k- kf) {k + li + li> + 1)}^/^ X 

x{k)-^''^{li C(^) C^'^) If) N''-^ {njljnili,nflfnidi>) , 



(8) 



(riiXinjXj Wh^soo^^"'"^'^'^]] ni>Ximj'Xj>) = -2 ■ 2"^ {2k + 1)^/^ [li ||c(*^)|| Z^/) x 
X {{k{k-hl)y^/'^{li{li-hl)-k{k-hl)-li' {k' + l)) X 

X + 1) N''~^ {njljnili,nj'ljini'lii) — kN^ {nilinjlj,ni'liinj'lji)^ — 

-2 {k {k + 1))^/^ 1/^"^ {nilinjlj,ni>li>njdj>)^ , 



(9) 



(niXiUjXj H^k+iki,na2i) m^Xim^Xf) = 2-2''^ {{2k + 1) {2k + 3) x 

X {k + li> -k){k-li + li> + l){k + li- li> + 1) {k -\- k -\- li> + 2)}^/^ X 

X (A: + 1)-^/^ (li ||C(*^)|| /i/) (ij If) N'' {nilinjlj,ni>li>nf{ 

The radial integrals in (8)-(10) are (see Glass and Hibbert [7]): 



'j') 



(10) 
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d 



Pi (ri) Pj (r2) ^r2— (n) Pf (ra) dndra- 



(11) 



Now we have all we need (the operators for tensorial structure and their submatrix elements) 
for obtaining the value of a matrix element of these operators for any number of open shells in bra 
and ket functions. This lets us exploit all advantages of the approach by Gaigalas et al [3]. 

The spin-spin and spin-other-orbit operators itself generally contain tensorial structure of several 
different types. Therefore the expression (1) must be used separately for each possible tensorial 
structure for performing spin-angular integrations according to [3]. Each type of tensorial structure 
is associated with a different type of recoupling matrix P (^Aj, Aj, A^-, A''™, A*^*^*, and with 

different matrix elements of standard tensorial quantities T (riiXi, njXj,n'^X[, n^A^-, A^*"", A'^^*, H, t) . 
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4 Conclusions 

The tensorial forms of the general secondly quantized spin-spin interaction operator (2) and spin- 
other-orbit interaction operator (7) and its submatrix elements ( for spin-spin interaction expres- 
sions (3), (4) and for spin-other-orbit expressions (8), (9) and (10)) are presented. In calculating 
its matrix elements between functions with u open shells this allows to exploit all the advantages 
of method by Gaigalas et al [3]: 

1. to obtain both diagonal and off-diagonalelements with respect to the configuration matrix 
elements in a unified approach, 

2. to use in practical applications the tables of submatrix elements of standard quantities, which 
here are both the coordinate representation and the occupation number representation ten- 
sorial operators, 

3. to apply the quasispin formalism the for occupation numbers parts and make use of it, 

4. to make use of having recoupling matrices simpler than in other known approaches. 
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